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Executive  Summary 


This  document  describes  a  specular  multipath  simulation  for  antennas  located  close 
to  a  sea  surface.  The  antennas  are  circularly  polarized  but  the  ray-trace  implemen¬ 
tation  readily  accommodates  any  other  polarizations.  The  transmitted  signal  St(£) 
arrives  at  the  receiver’s  antenna  along  the  direct  path  and  from  multiple  reflections 
off  the  sea  surface.  The  surface-level  antenna  turns  the  multipath  into  multiplicative 
noise — assuming  the  receiver  is  narrow  band.  Under  this  narrow-band  assumption, 
the  received  signal  Sn(t)  is  approximated  by  modulating  the  transmitted  signal  St(^) 
with  this  multiplicative  noise: 


SR(t )  =  a(t )  x  sT(t). 

The  multiplicative  noise  (a(t)}  is  determined  by  the  sea  surface,  the  elevation  angle 
of  the  transmitted  signal  path,  and  the  receiver’s  antenna.  By  sweeping  over  various 
realizations  of  the  sea  surfaces,  elevation  angles,  and  various  antenna  heights  and 
speeds,  the  systems  engineer  can  estimate  sea-surface  specular  multipath  effects  on 
the  surface-to-satellite  link. 
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Table  1:  Table  of  Notation. 


CP 

circularly  polarized 

RCP 

right-handed  circularly  polarized 

LCP 

left-handed  circularly  polarized 

Gt{0) 

gain  of  the  transmitter’s  antenna 

Gr(0) 

gain  of  the  receiver’s  antenna 

e 

elevation  angle 

9x 

elevation  angle  to  the  transmitter 

o9 

grazing  angle 

Oi 

angle  of  incidence 

9r 

angle  of  reflection 

r 

surface  reflection  coefficient 

r0 

plane  surface  reflection  coefficient 

To.y 

plane  surface  reflection  coefficient  for  vertical  polarization 

To  ,H 

plane  surface  reflection  coefficient  for  horizontal  polarization 

r0,c 

plane  surface  reflection  coefficient  for  circular  polarization 

To,x 

plane  surface  reflection  coefficient  for  cross-circular  polarization 

h-R 

height  of  the  receiving  antenna 

h+R 

upper  height  of  the  receiving  antenna  aperture 

h-R 

lower  height  of  the  receiving  antenna  aperture 

Xt 

wavelength  of  the  transmitted  wave 

pT 

polarization  of  the  transmitter’s  antenna 

Pr 

polarization  of  the  receiver’s  antenna 

E 

complex  electric  held 

Eh 

horizontal  component  of  E 

Ev 

vertical  component  of  E 

V 

polarization  efficiency 

sT(t) 

transmitted  signal 

SR{t ) 

received  signal 

C(r) 

height  of  the  sea  surface  at  r 

standard  deviation  of  the  sea  surface  height 

Hs 

significant  wave  height 

DSR 

Direct-to-Specular  Ratio 

QHA 

Quadrihlar  Helix  Antenna 

IX 


1  Time- Varying  Random  Channels 


Starting  in  the  early  1960s.  Bello  [7]  created  a  mathematical  framework  for  random, 
time-varying 1  filters  to  model  RF  channels.  This  time-varying  channel  model  maps 
the  transmitted  signal  sr(t)  to  the  received  signal  sn(t)  as  discussed  by  Bello  [7],  [8], 
[9],  [10]: 

/OO 

h(t,t-  r)sT(r)dr. 

-OO 

The  channel  model  h(t,  r)  is  called  the  input  delay-spread  function.  A  canonical 
example  is  the  Quadrature- Modulation  Fading  Simulator  [31],  [35,  pages  697-703]: 

N 

h(t,  r)  =  J2  an(t)S(r  -  Tn). 

n=  1 

In  this  model,  the  received  signal  is  the  sum  of  delayed  and  faded  versions  of  the 
transmitted  signal: 

N 

SR(t )  =  V  an(t)sT(t  -  Tn). 
n=  1 

Specihc  channels,  such  as  the  surface  scattering  channel,  are  described  by  the  de¬ 
lays  and  fading  processes.  In  particular,  if  the  delays  rn  are  not  resolvable  in  the 
bandwidth  fg  of  the  receiver,  say 


|t„|  <  10  Vs, 

then  the  channel  can  be  approximated  as  multiplicative  noise : 

SR(t )  «  |  V  X  ST(t).  (1) 

This  report  simulates  the  multiplicative  noise  caused  by  the  transmitted  signal  re¬ 
flecting  off  a  sea  surface  and  received  by  a  surface-level  antenna. 

Equation  1  is  obtained  by  making  several  simplifying  assumptions.  Section  2  in¬ 
troduces  the  model  assumptions  that  restrict  us  to  in-plane  scattering.  Consequently, 
we  only  will  scatter  off  a  one-dimensional  sea  surface  rather  than  a  two-dimensional 
sea  surface. 

Scattering  off  a  one- dimensional  surface  is  still  difficult — even  for  a  signal  bounc¬ 
ing  off  a  flat  surface,  there  is  a  mix  of  specular  and  diffuse  scattering.  Section  3 
illustrates  both  types  of  scattering  and  discusses  the  restriction  to  specular  scatter¬ 
ing.  The  restriction  to  specular  scattering  offers  a  simple,  qualitative  simulation  of 
the  multipath  noise. 

JKey  terms  are  italicized. 
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Section  4  reviews  several  quadrifilar  helical  antenna  patterns.  An  antenna  model 
is  developed  that  accounts  for  the  circular  polarization  on  the  direct  path  and  the 
cross-circular  polarization  that  arrives  from  the  sea-surface  reflections.  This  model 
can  be  coupled  to  the  flat-plane  scattering  to  estimate  the  direct-to- specular  ratio 
(DSR)  that  will  benchmark  the  simulations. 

Section  5  reviews  the  sea-surface  models.  The  models  are  two-dimensional  Gaus¬ 
sian  surfaces  obtained  by  filtering  the  sea-surface  spectrum.  Looking  into  the  wind 
of  a  fully  developed  sea  surface,  we  can  see  that  in-plane  scattering  is  a  reasonable 
approximation.  Therefore,  by  picking  off  a  “slice”  of  a  sea  surface  that  points  into 
the  wind,  we  have  a  one-dimensional  surface  appropriate  for  in-plane  scattering. 

Section  6  adapts  the  flat-plane  scattering  of  Section  3  to  local  tangent  planes  of  the 
one-dimensional  sea  surface.  That  is,  the  scattering  mechanism  are  little  flat  plates 
tangent  to  the  sea  surface  that  reflect  a  ray  from  the  transmitter  into  the  antenna 
aperture.  This  specular  multipath  noise  is  the  multiplicative  noise  in  Equation  1. 

Section  7  applies  this  specular  tangent-plane  scattering  to  several  sea  surfaces  and 
antennas.  These  simulations  show  that  the  specular  multipath  reflecting  off  a  smooth 
sea  surface  follows  the  two-path  DSR.  This  DSR  depends,  in  part,  on  the  receiving 
antenna.  For  the  low-angle  antennas  used  in  the  simulations,  the  bulk  of  the  specular 
multipath  arrives  from  the  lower  elevation  angles  (~  10°).  However,  as  the  sea  surface 
becomes  less  smooth,  the  specular  multipath  migrates  into  the  upper  elevation  angles. 

Section  8  emphasizes  that  this  specular  multipath  model  only  obtains  a  qualita¬ 
tive  model  for  Equation  1.  The  approximations  made  to  obtain  this  qualitative  model 
are  non-trivial.  Consequently,  the  temptation  for  endless  tweaking  is  unbounded. 
Tweaking  is  kept  in  check  by  knowing  how  the  approximations  intertwine — a  single 
small  “improvement”  does  not  “improve”  the  multipath  model.  Rather,  all  the  ap¬ 
proximations  must  be  simultaneously  tightened  before  the  multipath  model  registers 
significant  improvements.  Such  a  shift  in  the  modeling  fidelity  migrates  this  qualita¬ 
tive  model  into  a  full-up  scattering  model.  Consequently,  there  is  no  middle  ground 
in  scattering — either  a  simple  model  is  used  with  the  understanding  that  the  results 
are  qualitative  or  a  full-blown  propagation  model  is  used  that  accurately  captures 
both  the  local  scattering  off  the  sea  surface  and  the  long-range  propagation  through 
the  atmosphere  and  surface  ducts. 


2  In-Plane  Scattering 

Figure  1  illustrates  an  electric  field  scattering  off  a  two-dimensional  surface.  The  in¬ 
cident  electric  field  is  decomposed  into  a  vertical  component  determined  by  the  inci¬ 
dent  plane  and  the  horizontal  component  orthogonal  to  the  incident  plane  (right-hand 
rule).  Likewise,  the  reflected  electric  field  is  decomposed  into  a  vertical  component  in 
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the  reflected  plane  and  the  horizontal  component  orthogonal  to  the  reflected  plane.  If 
the  transmitter  is  aimed  orthogonal  to  the  wave  direction,  the  bulk  of  the  scattering 
will  stay  near  the  incident  plane.  Equivalent  descriptions  of  this  “in-plane  scattering” 
are  as  follows: 

•  The  scattering  is  such  that  purely  horizontal  or  purely  vertical  polarization  is 
not  depolarized. 

•  There  is  only  longitudinal  scattering — not  lateral  scattering. 

•  The  normal  to  any  reflecting  plane  lies  in  the  plane  of  incidence. 


Figure  1:  Coordinates  for  incident  and  reflected  polarizations  of  the  electric  field. 

In-plane  scattering  allows  the  use  of  a  one-dimensional  “slice”  of  a  sea  surface  for 
this  Phase  1  effort.  Phase  2  will  tackle  the  full  two-dimensional  sea  surface.  For 
one-dimensional  scattering,  assume  the  following: 

RF-1  Homogenous  atmosphere. 

RF-2  Distant  transmitter. 

RF-3  Surface- level  receiver. 

RF-4  No  surface-level  diffraction  or  aerosol  effects. 

RF-5  The  sea  surface  is  constant  in  each  plane  orthogonal  to  the  signal  direction. 

A  homogenous  atmosphere  has  a  constant  index  of  refraction  or,  equivalently,  a  con¬ 
stant  speed  of  propagation.  In  a  homogenous  atmosphere,  rays  travel  in  straight  lines. 
The  distant  transmitter  justifies  the  plane- wave  approximation  of  the  transmitter’s 
EM  field.  The  surface-level  antenna  justifies  a  flat-earth  approximation  for  ray  trac¬ 
ing.  The  absence  of  diffraction  over  the  waves  allows  ray  tracing  after  reflection. 
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Finally,  RF-5  forces  in-plane  scattering.  There  is  no  lateral  scattering.  Consequently, 
the  scattering  takes  place  in  a  two-dimensional  plane  reflecting  off  a  one-dimensional 
surface.  The  simplest  scattering  is  flat-surface  scattering. 


3  Flat-Surface  Scattering 

Figure  2  illustrates  the  ray  bundles  for  flat-surface  scattering.  A  transmitter  is  broad¬ 
casting  a  narrow-band  signal  sy(t)  at  wavelength  A t  and  polarization  pp. 2  The  distant 
transmitter’s  electric  held  is  approximated  by  plane  waves  traveling  in  the  direction 
shown  by  the  rays  with  elevation  angle  Op ■  The  receiver  has  an  antenna  located  at 
height  Hr. 


Figure  2:  Direct  path,  specular  reflection,  and  diffuse  scatter. 

The  direct  path  follows  free-space  propagation  to  the  antenna  to  arrive  delayed  in 
time  by  r\  and  scaled  by  the  complex  gain  g\  of  the  antenna  in  the  direct  path: 

SR.{t )  =  gisT(r  -  n). 

Because  the  signal  Spit)  is  narrow  band,  the  transfer  function  of  the  antenna  only 
multiplies  the  direct-path  signal  with  a  complex  gain.2 3 

The  plane  wave  also  illuminates  the  entire  surface.  Each  piece  of  the  surface 
projects  some  energy  into  the  antenna.  A  flat  surface  reflects  most  of  the  energy  in 
the  specular  direction  defined  by  Bello  [10]: 

The  direction  that  the  reflected  plane  wave  would  take  if  the  incident 
plane  were  reflected  from  a  mirror  located  at  the  mean  surface  height. 

2Pt  is  a  complex  number  characterizing  the  polarization  and  is  the  stereographic  projection  of 
the  Poincare  sphere  [30]. 

3There  is  no  loss  of  generality  for  ignoring  the  free-space  loss  and  scaling  only  by  the  antenna. 
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The  specular  reflection  is  modeled  as  a  ray  that  travels  to  the  antenna.  However, 
the  surface’s  reflection  coefficient  F  reduces  power,  changes  phase,  and  modifies  the 
polarization.  For  a  flat,  mirror-like  surface,  the  received  signal  is  the  sum  of  the  direct 
path  and  specular  reflection: 

SR(t)  =  gisT{r  -  ti)  +  g2FsT{T  -  r2), 

where  g 2  is  the  complex  gain  of  the  antenna  in  the  specular  direction.  Consequently, 
the  input  delay-spread  function  h(t,r)  has  the  form  determined  by  Matthews  [28,  Eq. 
2.25],  and  Parsons  [36,  Eq.  2.19]: 

h(t,r)  =  gi6(r  -  n)  +  g2V 6(r  -  r2)  .  (2) 

direct  path  specular  reflection 

In  this  flat-surface  model,  the  signal  does  not  change  if  the  surface  is  moved  horizon¬ 
tally.  This  model  can  be  extended  to  slightly  rough  surfaces  by  incorporating  diffuse 
scatter.  For  example,  Figure  3  shows  a  rough  surface  breaking  the  specular  reflec¬ 
tion  into  non-specular  directions  or  into  diffuse  scatter ,4  The  Fraunhofer  criterion  or 
Rayleigh  criterion  quantifies  the  surface  roughness  in  terms  of  the  standard  deviation 
oy-  of  the  sea  surface. 


Figure  3:  Specular  reflection  and  diffuse  scatter  off  a  slightly  rough  surface. 

The  surface  is  considered  smooth  and  produces  a  strong  specular  reflection  pro¬ 
vided  [15,  Eq.  1],  [36,  Eq.  2.29]: 

The  RMS  phase  difference  between  two  rays  re- 
^  J_  x  ^ T  fleeted  at  two  different  heights  on  the  surface  must 

^  32  sill(6fp)  be  smaller  than  tt/8  in  the  far  field  to  coherently" 

combine. 

4John  Rockway,  Akira  Ishimaru,  Yasuo  Kuga,  Seung-Woo  Lee  [2005]  Statistical  Green’s  Functions 
for  Applications  to  Rough  Surface  and  Random  Media  Scattering,  Object  Interaction  and  Detection, 
University  of  Washington  Electromagnetics  and  Remote  Sensing  Laboratory,  Power  Point  Presenta¬ 
tion. 
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Figure  4  displays  the  smooth-surface  region  (the  blue  region)  as  a  function  of  angle 
and  surface  deviation.  Because  the  sea  surfaces  in  this  simulation  typically  have 
cr^  <  0.5  meters,  the  figure  shows  that  the  flat-surface  model  starts  to  break  down 
when  the  transmitter’s  elevation  angle  exceeds  5  degrees.  Moreover,  the  Rayleigh 
criterion  does  not  account  for  the  blocking  of  the  rays  that  occurs  at  the  low  grazing 
angles. 


grazing  angle  (deg) 


Figure  4:  Rayleigh  criterion  (blue)  at  1-meter  wavelength  in  the  Ot-vc,  plane. 


A  typical  model  uses  a  Gaussian  random  process  {g{t)}  to  approximate  diffuse 
scattering  [26],  [11],  [29]: 

h(t,r)  =  9iS(r  -  ri)  +  ^2r^(r-r2)  +  g(t)6{r  -  r3)  .  (3) 

' - „ - '  - - . - '  - - . - ' 

direct  path  specular  reflection  diffuse  scatter 

The  following  sections  unpack  the  time  delays,  reflection  angles,  antenna  gains,  and 
reflection  coefficients  for  this  model.  The  final  section  specializes  this  model  to  spec¬ 
ular  scattering  into  a  surface-level  antenna. 
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3.1  Time  Delays 


For  the  receiver  simulation,  there  is  no  loss  of  generality  to  work  with  the  relative 
delays  by  taking  T\  —  (1.  In  Figure  5,  provided  the  transmitter’s  elevation  angle 


Figure  5:  Time-Delay  geometry  for  elevation  angle  9R  <  45°. 

9t  <  45°,  the  delay  between  the  direct  path  and  reflected  path  is 

r2  =  c~l(dR  -  dT)  =  2hRsm(9T).  (4) 

When  the  transmitter’s  elevation  angle  9R  exceeds  45  degrees,  the  geometry  is  differ¬ 
ent,  but  the  delay  formula  is  the  same.  Recalling  that  1  ns  ~  1  foot,  a  3-meter  mast 
exhibits  a  maximum  delay  of  approximately  18  ns.  Figure  6  shows  that  the  small 
elevation  angles  have  time  delays  of  approximately  1  nanosecond.  A  receiver  with  a 
bandwidth  of  100  kHz  cannot  resolve  these  delays  for  signal  enhancement.  Although 
the  simulations  show  that  longer  delays  are  possible  because  of  distant  reflections  off 
the  sea  surface,  even  these  longer  delays  do  not  exceed  a  few  hundred  nanoseconds. 
The  narrow-band  receiver  still  cannot  resolve  these  longer  paths.  Consequently,  the 
zero-delay  of  Equation  1  still  applies. 

3.2  Antenna  Gains 

If  the  polarization  pR  of  the  transmitted  signal  matches  the  polarization  pR  of  the 
receiving  antenna,  the  magnitude  of  the  gains  on  the  direct  and  reflected  paths  are 

\gi\  =  Gr(9t )^2,  |  <72 1  =  Gr{— 9t)1^2, 

where  the  receiving  antenna’s  gain  GR(9)  is  a  function  of  elevation  angle. 

To  handle  the  polarization  mismatch  between  the  receiving  antenna  with  polariza¬ 
tion  pR  and  a  single  transmitted  wave  with  polarization  py,  the  polarization  efficiency 
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Figure  6:  Time  delays  for  a  short  mast  and  flat  reflections. 


7j  modifies  the  gains  as  determined  by  Beckmann  [6,  page  186]: 

ltd  I  =  Gr^t)1^2  x  v(Pr(Qt),PtY ^  ■ 

Here,  p(pr,Pt)  is  the  ratio  of  the  power  delivered  to  the  antenna  from  a  wave  with 
polarization  px  to  the  maximum  power  delivered  to  the  antenna  from  a  wave  with 
polarization  pR  and  computed  as 

,  Y+PrPt\ 2 

'n{pR,PT)  (1  +  |ph|2)(1  +  |pt|2)' 

When  multiple  waves  of  various  polarizations  arrive  at  the  antenna,  the  antenna 
functions  as  an  electrical  transducer  mapping  the  received  waves  into  voltages  at  the 
antenna’s  terminal.  The  antenna’s  complex  effective  length  was  conceived  by  Sinclair 
[40]  to  model  this  transducer  and  effectively  complexifies  the  gains.  When  only  the 
antenna  gain  is  available,  this  limitation  is  acknowledged  by  using  the  gain  Gr{9). 

3.3  Reflected  Path 

The  ray  that  coherently  reflects  off  the  sea  surface  loses  power  and  is  depolarized. 
The  power  loss  and  depolarization  are  encoded  in  the  factorization  of  the  reflection 
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coefficient  [10],  [33,  pages  19-22],  [25,  pages  95-125]: 

r  =  Dr0n,  (5) 

where  D  denotes  the  divergence,  7 Z  denotes  the  roughness  coefficient,  and  r0  is  the 
plane  surface  reflection  coefficient. 

Divergence:  The  divergence  0  <  D  <  1  [25,  page  98]  “is  equivalent  to  a  purely 
geometric  factor  that  describes  additional  spreading  of  a  beam  of  rays  due  to  reflection 
from  a  spherical  surface.” 

Roughness:  As  the  sea  surface  changes  from  a  glassy  “mirror-like”  surface  to  a  wind- 
ruffled  chop,  the  specular  reflection  of  any  ray  decreases.  This  decrease  is  modeled  by 
the  surface  roughness  coefficient  1Z  derived  from  experimental  observations  [4],  [5], 
[20],  [25,  pages  122-125],  [13,  Eq.  2.62]: 

K(0g)  =  e'g2/2/0(g2/2);  g  =  47rsin(0s)^,  (6) 

where  g  is  the  “electrical  roughness  parameter”  [10]  or  “apparent  ocean  roughness” 
[4],  the  grazing  angle  6g  is  the  angle  between  the  ray  and  the  local  tangent  plane 
( 9g  =  9t  in  the  plot),  oy  denotes  the  standard  deviation  of  the  height  fluctuation  of  the 
surface  [33,  page  59],  [25,  pages  50-54],  and  A t  is  the  wavelength  of  the  transmitted 
waveform. 


Plane  Surface  Reflection  Coefficients:  The  reflection  coefficient  To  for  a  plane  sur¬ 
face  depends  on  the  complex  permittivity  e  of  the  terrain  and  the  polarization.  For 
sea- water  at  20  degrees  Centigrade  and  frequencies  less  than  1  GHz,  the  complex 
permittivity  is  [13,  Eq.  3.1]: 

e  =  e'  —  j60a\T  =  70  —  j'300At- 


The  reflection  coefficient  for  vertical  polarization  [25,  page  99],  [13,  page  53]: 

esin(0s)  —  Je  —  cos  (9g)2 

To,v  =  - ; - 7 

esin  (9g)  +  Je  —  cos  (9g)2 


(7) 


where  the  grazing  angle  9g  is  the  angle  between  the  ray  and  the  local  tangent  plane 
(9g  =  9t  in  the  plot).  For  horizontal  polarization,  the  reflection  coefficient  is  [21,  page 
396], 

^  _  sin (dg)  -  \Je  —  cos(9g)2  ^ 

t  0,H  ~  - ,  ■ 

sin (9g)  +  Je  -  cos (9g)2 

Figures  7  and  8  plot  the  reflection  coefficient  for  seawater  as  a  function  of  the  grazing 
angle.  The  scales  are  different  for  clarity  and  show  that  the  horizontal  reflection 
coefficient  for  seawater  is  essentially  —1. 
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Figure  7:  Vertical  reflection  coefficient  for  seawater. 


Reflection  Coefficient:  Horizontial  Polarization;  SeaWater  @  20  C 


Figure  8:  Horizontal  reflection  coefficient  for  seawater. 
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The  problem  is  to  use  these  vertical  and  horizontal  reflection  coefficients  to  com¬ 
pute  the  reflection  coefficient  for  any  elliptical  polarization.  Of  the  several  approaches, 
a  beautiful  “change-of-basis”  is  found  in  Stutzmann  [42],  The  real- valued  instanta¬ 
neous  electric  field  E ft)  can  be  written  as 

m  =  COs(uff)£AUff  +  cos(o )t  +  8)Ey  Uy,  (9) 


where  Uy  and  u#  are  the  orthogonal  unit  vectors  in  the  vertical  and  horizontal 
directions.  A  left-handed  circularly  polarized  (LCP)  wave  has  Ey  =  Eh  and  <5  =  +90° 
[42,  Eq.  2.30]  so  that 


E  Lit) 


El 

a/2 


{cos(aff)u# 


sin(off)uy}  . 


A  right-handed  circularly  polarized  (RCP)  wave  has  Ey  =  Eh  and  S  =  —90°  [42,  Eq. 
2.31]  so  that 

E 

E^(t)  =  — f=  {cos(ut  +  8')uh  +  sin  (off  +  5')uy}  , 

V2 

where  8'  models  a  phase  difference  between  El{8)  and  Enft).  Any  electric  field  Eft) 
can  be  represented  in  the  Eh  ft)  and  Ey(£)  basis.  Likewise,  any  electric  field  can  be 
represented  in  the  El  ft)  and  E  nft)  basis.  The  complex  electric  field  provides  a  conve¬ 
nient  change-of-basis  from  linear  to  circular  polarization.  The  instantaneous  electric 
field  Eft)  links  to  the  phasor  electric  field5  intensity  E  as  provided  by  Stutzmann  [42, 
Eq.  3.11]: 

Eft)  =  $t[Eejujt\. 

The  phasor  electric  field  may  be  decomposed  as  provided  by  Stutzmann  [42,  Eq.  3.12] 


E  —  EhVLh  +  Ey  Uy , 


where  Eh  and  Ey  are  now  their  complexified  versions  in  Equation  9.  The  complexified 
components  of  circular  polarization  relate  to  the  complexified  components  of  the  linear 
polarizations  as  provided  by  Stutzmann  [42,  Eq.  8.31]: 


'  Ey  ' 

1 

'  i 

1  ' 

'  El  ' 

.  Eh  _ 

“72 

_  3 

~j  _ 

_Er 

'  El  ' 

1 

'  1 

~j 

'  Ey  ' 

.  er  . 

“71 

1 

j  _ 

Eh 

Assumption  RF-5  lets  us  claim  in-plane  scattering  or  that  the  scattering  is  such  that 
purely  horizontal  or  purely  vertical  polarization  is  not  depolarized.  That  is,  the  in¬ 
cident  electric  field  E*  is  mapped  to  the  reflected  electric  field  Er  as  provided  by 
Stutzmann  [42,  Eq.  8.29]: 


Ery 

To,v  0 

Ej.y 

Er,H 

o 

o 

Ei,H 

5The  Fourier  transform  of  ER). 
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The  change-of-basis  converts  this  scattering  of  the  linear  polarizations  to  circular 
scattering  provided  by  Stutzmann  [42,  Eq.  8.34]: 


1 

i 

o 

o 

To.x 

Ei,L 

Er,R 

r0,x 

To,c 

EhR 

where  the  circular  and  cross-circular  reflection  coefficients  are  [42,  Eq.  8.35] 

To.c  =  2  {IVv  —  r 0,ij}  , 

To.x  =  2  (ro,v  +  r0  ,h}  , 

and  plotted  in  Figures  9  and  10,  respectively. 


Reflection  Coefficient:  Circular  Polarization:  SeaWater  @  20  C 


grazing  angle  (deg) 


Figure  9:  Circular  reflection  coefficient  for  seawater. 
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Reflection  Coefficient:  Cross  Polarization;  Sea  Water  @  20  C 


grazing  angle  (deg) 


Figure  10:  Cross-circular  reflection  coefficient  for  seawater. 


3.4  Diffuse  Scatter 

As  the  sea  surface  changes  from  a  glassy  surface  to  a  wind-ruffled  chop,  the  specular 
reflection  decreases  and  the  diffuse  scatter  increases.  The  diffuse  scatter  is  typically 
modeled  by  a  zero-mean,  narrow-band,  complex-valned  Gaussian  random  process 
(g(t)}.  Consequently,  only  the  spectral  shape  and  variance  a2g  need  be  specified. 
Beard  [5],  [4]  obtains  a  heuristic  model  for  a2  using  the  surface  roughness,  whereas 
Ryan6  argues  from  conservation  laws  that 

a*  =  Gf\r0\\i-n2).  (io) 

Here  Gf  denotes  a  “forward-scatter”  gain.  Regardless  of  the  approach,  a2  is  a  func¬ 
tion  of  the  surface  roughness  TZ — assuming  the  model  implicit  in  Equation  3.  The 
assumption  for  this  model  are  discussed  in  the  next  section. 


6Personal  communication  from  Frank  Ryan,  SSC  San  Diego. 
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3.5  Specular  Scattering  at  a  Surface-Level  Antenna 

A  qualitative  model  of  specular  scattering  into  a  surface-level  antenna  is  constructed 
from  the  preceding  computations.  This  model  captures  those  effects  of  specular  scat¬ 
tering  that  generalize  to  tangent-plane  scattering  off  a  sea  surface. 

Placing  the  antenna  near  the  a  relatively  smooth  surface  allows  us  to  ignore  the 
time  delays,  the  divergence,  the  surface  roughness  coefficient,  and  the  diffuse  scatter. 
However,  placing  the  antenna  within  a  wavelength  of  the  surface  distorts  the  antenna 
pattern. 

Regarding  the  time  delays,  Section  3.1  shows  that  the  low  antenna  induced  time 
delays  by  a  few  nanoseconds.  The  sea-surface  simulations  show  that  even  the  distant 
reflections  delay  the  signal  by  only  a  few  hundred  nanoseconds.  If  the  receiver  has 
a  bandwidth  of  only  100  kHz,  these  nanosecond  delays  cannot  be  resolved.  Conse¬ 
quently,  the  time  delays  are  ignored. 

Regarding  the  divergence  D ,  the  sea-surface  simulations  show  reflections  are  lim¬ 
ited  to  a  few  kilometers.  Consequently,  the  flat  earth  is  an  excellent  approximation 
and  D  =  1  [21,  page  113]. 

Figure  11  plots  the  surface  roughness  coefficient  TZ  of  Equation  6  as  a  function 
of  the  standard  deviation  oy  of  the  surface.  The  wavelength  is  At  =  1  meter.  Each 
curve  is  plotted  for  a  fixed  elevation  angle:  Ox  =  5°,  10°,  30°.  Because  the  sea  surface 
simulations  have  oy  <  0.5  meters,  the  surface  roughness  coefficient  7 Z  approximately  1 
for  elevation  angles  6t  <  10°.  Consequently,  its  effect  on  surface  reflection  coefficient 
T  of  Equation  5  is  ignored. 

For  the  larger  angles  and  the  rougher  seas,  7 Z  is  fit  to  widely  varying  experimental 
data.  From  Long  [25,  pages  122-124]: 

•  Reflection  data  are  often  given  with  little  precision. 

•  Evaluation  of  T  is  only  an  estimate. 

•  Values  for  oy  are  at  best  guesses. 

Moreover,  the  form  of  7 Z  neglects  edge  and  shadowing  effects.  Consequently,  we  take 
TZ  =  1  up  to  30°,  noting  that  our  sea-surface  scattering  will  include  shadowing  and 
the  roughness  caused  by  multiple  reflectors  within  the  antenna’s  footprint. 

The  following  observations  justify  omitting  the  diffuse  scattering  term  of  Equa¬ 
tion  3.  First,  the  diffuse  scatter  {g(t)}  approximates  the  effect  multiple  reflectors 
within  the  antenna’s  footprint.  Because  the  model  we  develop  includes  these  multi¬ 
ple  reflectors,  some  diffuse  scattering  is  actually  included.  Second,  Equation  10  shows 
the  diffuse  scatter  is  small  when  TZ  ~  1.  The  surface  roughness  TZ  depends  on  the 
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Surface  Roughness  Coefficient  [Eq.  2.62,  Boithias] 


Figure  11:  Surface  roughness  coefficient. 


variance  a2  of  the  sea  surface  £(r).  This  variance  depends  on  the  wind  speed  U  from 
Equations  13  and  14.  Figure  12  plots  the  diffuse  scatter  power  a2  at  300  MHz  as 
a  function  of  wind  speed  U  for  grazing  angles  6 t  =  5°,  10°,  and  30°.  The  forward 
scatter  gain  is  Gf  =  1/2.  The  plot  shows  that  omitting  a2  induces  errors  10%  or  less, 
provided  that  the  wind  speed  is  less  than  10  m/s  or  that  the  grazing  angle  is  less  than 
10  degrees. 

With  negligible  time  delays,  surface  roughness  coefficient  1Z  =  1,  and  diffuse 
scatter  ignored,  the  input  delay-spread  function  for  this  Flat-Surface  Scattering  is 
well-approximated  by  the  direct  path  and  a  single  reflection: 

h(t,  r)  =  {Gr(Ot)1^2  +  Gr(— ^t)1^2To,c(^t)}  x  (11) 

Section  6  applies  this  model  to  develop  tangent  plane  scattering  used  for  the  multipath 
simulations. 

In  summary,  this  Flat-Surface  Scattering  Model  of  Equation  11  omits  the  follow¬ 
ing: 


•  Delays:  Nanosecond  delays  are  not  resolvable  with  a  kHz  bandwidth. 

•  Diffuse  scatter  (7 Z  =  1):  The  flat  surface  is  assumed  to  be  smooth. 

•  Curved  Earth  ( D  =  1):  The  surface-level  antenna  allows  us  to  assume  a  flat 
earth  [21,  page  113]. 

•  Scintillation:  No  fast  fading  on  the  rays. 
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Figure  12:  Diffuse  scatter  power  a2g  as  a  function  of  wind  speed. 

•  Shadowing  and  blocking:  Not  relevant  because  of  flat-surface  assumption. 

•  Loss:  Propagation  loss  on  the  reflected  path  is  ignored. 

•  Antenna  mismatch:  The  transmitting  and  receiving  antennas  have  the  same 
circular  polarization. 

The  next  section  introduces  the  antenna  patterns  Gr(9 )  and  extends  Equation  11  to 
account  for  the  antenna  mismatch. 
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4  Antenna  Patterns 


The  simulations  use  a  quadrifilar  helix  antenna  as  the  receiving  antenna.  Figure  13 
is  a  picture  of  a  typical  quadrifilar  helix  antenna.  Section  4.1  displays  a  text-book 
pattern  produced  by  a  quadrifilar  helix  antenna.  This  pattern  has  a  simple  cardiod 
shape  so  is  an  excellent  baseline  case.  Section  4.2  reports  on  measured  patterns  for  a 
right  circularly  polarized  (RCP)  quadrifilar  helix  antenna.  These  patterns  are  more 
complex  compared  to  the  model  quadrifilar  helix  antenna.  The  plots  also  show  the 
antenna’s  response  to  both  right  and  left  circular  polarizations.  Section  4.4  makes 
the  final  adjustments  in  the  multipath  model  to  include  both  the  RCP  and  LCP 
measurements. 


Figure  13:  Electronically  Reconhgurable  SATCOM  Antenna  (courtesy  of  Toyon  Re¬ 
search  Corporation). 
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4.1  Baseline  Quadrifilar  Pattern 


Kilgus  [22],  [23]  analytically  investigated  the  quadrifilar  helix  and  determined  that 
this  antenna  exhibits  a  cardiod  pattern.  Figure  14  shows  such  a  pattern  as  measured 
by  Sainati  [38].  Because  the  sea-surface  reflections  arrive  in  the  lower  hemisphere, 
the  lower  hemisphere  determines  the  multipath  effect.  For  this  reason,  the  cardiod 
pattern  is  the  baseline  case. 


Quadrifilar  Helical  Antenna  Pattern  at  250  MHz 


0 

OdB 


Figure  14:  Measured  gain  as  a  function  of  elevation  angle  for  a  resonant  quadrifilar 
helix  [38]. 
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4.2  Electronically  Reconfigurable  Antenna  (ERA)  Patterns 

Figures  15  and  16  show  the  Electronically  Reconfigurable  Antenna  (ERA)  patterns. 
The  antenna’s  RCP  response  is  plotted  in  blue.  Plotted  in  red  is  the  response  to  LCP. 
The  low-beam  ERA  peaks  at  approximately  10  degrees  elevation.  The  patterns  are 
unexpectedly  asymmetric.  The  asymmetries  may  be  caused  by  coaxial  cable  inter¬ 
ference,  the  ground  plane  being  too  small,  and  off-center  positioning.  In  the  upper 
hemisphere,  the  RCP  gain  exceeds  the  LCP  gain  by  10  dB.  In  the  lower  hemisphere, 
the  LCP  exceeds  the  RCP.  The  reflections  from  the  sea  surface  arrive  in  the  lower 
hemisphere  and,  at  low  grazing  angles,  the  sea  surface  changes  the  polarization  from 
RCP  to  almost  LCP.  Consequently,  the  LCP  measurement  is  needed  to  register  the 
LCP  multipath.  Section  4.4  shows  how  to  use  both  patterns  to  model  the  multipath. 


iUHF  Helix:  low  beam  1  @  250  MHz;  measured  12/05/00 
0° 

0  dB 


angle  from  zenith  in  athwart  plane 


Figure  15:  Toyon’s  low-beam  ERA  at  250  MHz  (courtesy  of  Toyon  Research  Corpo¬ 
ration). 
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iUHF  Helix:  low  beam  1  @310  MHz;  measured  12/05/00 
0° 

OdB 


angle  from  zenith  in  athwart  plane 


Figure  16:  Toyon’s  low-beam  ERA  at  310  MHz  (courtesy  of  Toyon  Research  Corpo¬ 
ration). 
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4.3  OE-538  Patterns 


Measurements  for  the  OE-538  antenna  were  furnished  by  the  Naval  Undersea  Warfare 
Center  (NUWC).  The  patterns  are  plotted  in  Figures  17,  18,  19,  and  20.  The  patterns 
are  azimuthal  (</>)  slices  of  the  electric  held  magnitudes  \Ec/)\  and  \Eg\  plotted  against 
the  elevation  angle  6.  These  E-held  gains  are  measured  at  the  base  connector  of  the 
OE-538  Antenna  Assembly.  Consequently,  the  gains  are  reduced  by  1.55  dB  because 
of  losses  in  the  internal  cabling,  diplexers,  and  switches7. 


OE-538  High  Angle  IE  I  @  300  MHz 

0  5  dB 


Figure  17:  E0  at  300  MHz  for  High- Angle  OE-538  antenna  (courtesy  of  NUWC). 


7Porco,  Pietro  G.  [2004]  OE-538  SATCOM  Antennae  Gain  Footnotes,  NUWCDIVNPT. 
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OE-538  High  Angle:  \E  \  @  300  MHz 

0 


0° 


Figure  18:  Eg  at  300  MHz  for  High- Angle  OE-538  antenna  (courtesy  of  NUWC). 
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OE-538  Low  Angle:  \E\  @  300  MHz 

0  5  dB 


Figure  19:  at  300  MHz  for  Low- Angle  OE-538  antenna  (courtesy  of  NUWC). 
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OE-538  Low  Angle:  \E  \  @  300  MHz 

0 


0° 


Figure  20:  Eq  at  300  MHz  for  Low-Angle  OE-538  antenna  (courtesy  of  NUWC). 
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The  RCP  and  LCP  responses  can  only  be  estimated  when  the  phase  is  missing. 
The  phasor  electric  held  may  be  decomposed  as  [42,  Eq.  3.12] 

E  =  EffUff  +  Eylly. 

We  estimate  the  horizontal  and  vertical  components  by  assuming  RCP  or  6  =  —90°: 

Eh  =  \E<j,\ 

Ev  =  \Ee\e~j^2 

Convert  to  CP  fields  as  provided  by  Stutzman  [42,  Eq.  8.31]: 


'  El  ' 

1 

'  i 

-3 

'  Ey  ' 

.  er  . 

“71 

i 

3  . 

Eh 

Figure  21  plots  \E^\  and  \Eq\  at  (j)  =  0°  azimuth.  Characteristic  of  a  CP  antenna,  both 
fields  have  nearly  equal  magnitude.  Converting  these  fields  to  Er  and  El  produces 
the  plots  in  Figure  22.  This  plot  shows  the  estimated  RCP  in  blue  and  the  estimated 
LCP  response  in  red. 

OE-538  Low  Angle:  E  fields  @  300  MHz  and  0  azimuth 


Figure  21:  E  fields  at  300  MHz  for  Low- Angle  OE-538  (courtesy  of  NUWC). 
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OE-538  Low  Angle:  Estimated  ER  and  ER  @  300  MHz  and  0  azimuth 


Figure  22:  Estimated  RCP  and  LCP  gain  at  300  MHz  for  Low-Angle  OE-538. 
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4.4  Antenna  Mismatch 


When  a  RCP  ray  hits  a  flat  sea  surface  at  a  low  grazing  angle,  two  rays  are  reflected: 

•  A  weak  RCP  ray  (Figure  9) 

•  A  strong  LCP  ray  (Figure  10) 

Arriving  at  an  ideal  RCP  antenna,  the  LCP  rays  would  not  generate  any  voltage  in 
the  antenna.  However,  the  antenna  plots  show  that  a  real-world  RCP  antenna  is  not 
ideal — the  real-world  RCP  antenna  factors  into  an  ideal  RCP  antenna  and  an  ideal 
LCP  antenna.  Figure  23  illustrates  this  antenna  and  reflection  where  the  received 
signal  is  the  voltage  sum  from  both  ideal  antennas:  The  weak  RCP  ray  exciting 
the  ideal  RCP  antenna  and  the  stronger  LCP  ray  exciting  at  the  weaker  ideal  LCP 
antenna. 


Non-Ideal  Antenna 


Figure  23:  RCP  reflecting  from  sea  surface  and  received  by  non-ideal  RCP  antenna. 

Figure  23  extends  Equation  11  to  account  for  both  RCP  and  LCP  reflections  from  a 
single  ray: 


direct  path  specular  RCP  reflection 

h (t ,  7~)  =  (Gr.RCP^t)1/2  +  Grircp(  — 

+  Gr,LCp(_ ^T^^To.X^t)  }  X  <5(t).  (12) 

^ - V - ' 

specular  LCP  reflection 

However,  a  wind-roughened  sea  surface  offers  many  reflection  points.  Sea  surfaces  are 
introduced  in  Section  5.  Their  reflection  points  are  modeled  in  Section  6. 
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5  Sea  Surfaces 


The  sea-surface  waves  are  assumed  to  be  generated  by  a  wind  blowing  across  an  open 
fetch  of  ocean.  An  excellent  description  of  the  wind  action  on  the  sea  surface  and  the 
mechanisms  producing  the  wave  types  is  provided  by  the  Department  of  Geosciences, 
University  of  Houston  [14]: 

A  gentle  breeze  blowing  across  the  surface  causes  stresses  that  are  ex¬ 
pressed  in  capillary  waves.  These  are  tiny  ripples  with  rounded  crests, 
V-shaped  troughs,  and  wavelengths  less  than  1.74  cm.  They  die  down 
immediately  if  the  wind  stops  blowing,  by  the  restoring  force  of  capillar¬ 
ity  (surface  tension)  . . .  These  ripples  roughen  the  surface,  providing  more 
vertical  surface  for  the  wind  to  push.  As  the  continues  to  blow,  the  ripples 
grow  larger  and  are  tranformed  into  gravity  waves  . . .  The  wavelength  is 
generally  15  to  35  times  the  height.  If  the  wind  stops  blowing,  these  waves 
die  down  by  the  restoring  force  of  gravity  ...  As  winds  continue  to  blow, 
these  gravity  waves  grow  in  size.  Wave  height  increases  more  rapidly  than 
wavelength.  The  shape  changes,  with  crests  becoming  more  angular  and 
troughs  rounded.  Eventually  the  wave  reach  a  maximum  size  for  that  par¬ 
ticular  wind  velocity,  defined  as  the  point  at  which  the  energy  supplied 
by  the  wind  is  equaled  by  the  energy  lost  by  breaking  waves  (white-caps). 

This  is  called  the  fully  developed  sea. 

Figure  24  lays  out  the  sea-surface  coordinates  and  parameters.  The  wind  with 
mean  speed  U  blows  into  angle  6  over  a  sea  at  temperature  Tw  for  a  duration  measured 
by  the  inverse  wave  age  D.  As  indicated  in  Figure  24,  spacial  positions  are  elements 
of  R2  denoted  as 

Tl 

r  = 

r2 

The  height  of  the  sea  surface  is  denoted  as  C(rV)  and  is  a  stochastic  random  held 
that  is  a  function  of  position  r  and  time  t.  The  duration  of  the  signals  are  so  small 
that  the  sea  surface  may  be  considered  stationary  or  frozen  in  place.  Consequently, 
the  time  dependence  is  omitted.  Typically,  the  sea-surface  height  <C(r)  is  assumed  to 
be  wide-sense  stationary  with  a  mean  height  of  zero 

o  =  £[C(r)]  M, 

a  variance  of 

<A2  =  £[C(r)2]  [m2], 

and  covariance  that  depends  only  the  spatial  difference: 

R(dr)  =  £[C(r/  +  r)C(r')]  [m2]. 
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Figure  24:  Sea-surface  parameters. 

The  significant  wave  height  Hs  is  the  mean  of  the  highest  one-third  of  the  waves 
present  in  a  sea  [24,  page  302],  Massel  derives  [27,  pages  139-141] 


assuming  the  sea  is  narrow-band,  Gaussian,  and  fully  developed.  A  fully  developed 
sea  requires  both  time  and  area  to  deliver  wind  energy  to  raise  the  waves  [18]:  a  wind 
10  m/s  needs  to  blow  with  a  duration  of  18  hours  over  about  320  km  of  ocean  The 
significant  wave  height  Hs  links  to  the  sea-surface  parameters  derived  by  Massel  [27, 
page  81]: 

U2 

Hs  =  0.2092  —  [m],  (14) 

where  g  is  the  gravitational  constant.  The  mean  wind  speed  U  denotes  the  measure¬ 
ment  taken  at  10  meters  height  [17,  Eq.  17].  Comparing  Hs  from  Equation  14  to  the 
statistical  estimate  from  Equation  13  benchmarks  sea-surface  simulations. 

The  sea  surfaces  are  simulated  by  the  wavenumber  spectrum.  <h(k).  The  wavenum¬ 
ber  spectrum  <h(k)  is  the  Fourier  transform  of  the  covariance  of  the  sea  surface: 

/oo  roo 

/  e+i{fc^-^>$(A:1,A:2)dfc1dfc2, 

-OO  J  — OO 

where  the  integral  is  interpreted  in  Generalized  Random  Process  Theory  [19],  [1]. 
The  wavenumber  vector  is  denoted 
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Figure  25  shows  the  wavenumber  spectrum  down  to  50  dB.  The  spectrum  strongly 
peaked  at  the  peak  wavenumber  kp. 

<D(k)  [dB]:  U—\0  m/s,  e=0°,  7W=15  C,  Q=0.83 


Figure  25:  Wavenumber  spectrum  (3-D  display). 

Associated  with  the  peak  wavenumber  kp  is  its  peak  wave  that  travels  at  speed  cp. 
The  inverse  wave  age  Q  measures  how  closely  the  wind  is  aligned  with  the  dominant 
wave  [41]: 

>  0.8333  growing  waves  or  young  sea 
<  0.8333  fully  developed  sea  or  stable  swells  ’ 

where  A 6  is  the  angle  between  the  wind  and  dominant  wave  direction.  Given  the 
inverse  wave  age,  the  spectral  features  can  be  computed  for  both  a  young  sea  and  a 
fully  developed  sea.  For  example,  the  peak  wavenumber  and  significant  wave  height 
for  a  young  sea  are  modified  by  Zavoronty  and  Voronovich  [43]: 

a  U2 

h  =  D2-^-;  Hs  =  0.1540  hT1'65  — , 

U 2  g 

These  models  derive  from  the  dependence  of  wavenumber  spectrum  <h(k)  on  the 
inverse  waveage.  Table  2  lists  one  choice  of  variables  that  model  $(k).  Variables 
below  the  line  never  change  for  the  simulations.  Typically,  above  the  line,  only  the 
wind  speed  varies.  Figure  26  presents  the  spectrum  parameterized  by  the  nominal 
values  of  Table  2. 


Q  =  —  cos(A  9) 
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Table  2:  Gravity  and  capillary  spectral  variables  for  <f>(k). 


Variable 

Nominal  Values 

Description 

U 

10 

mean  wind  speed  at  10  meters  (m/s) 

^max 

0 

mean  wind  direction  (degrees) 

T 

L  w 

15 

seawater  temperature  (C) 

n 

0.83 

inverse  wave  age  (dimensionless) 

/slick 

1 

slick  capillary  wave  suppression  factor 

s 

35.0 

Salinity  (psu) 

Pw 

1024.8 

sea- water  density  (kg/m3) 

Ts 

0.074 

sea- water  surface  tension  (N/m) 

Pa 

1.225 

density  of  air  (kg/m3) 

9 

9.7976 

gravitational  constant  (m/s2) 

hi 

0.42 

Karman  constant 

<f>(  k)  [dB] 


-1  -0.5  0  0.5  1 

fcj( m"1):  U=  10  m/s;  9=0°;  T  =15  C;  0=0.83 


Figure  26:  Wavenumber  spectrum  (top  view). 
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Because  the  spectral  peak  is  difficult  to  see  in  the  linear  scale,  Figure  27  plots 
“slices”  of  the  spectrum  $(fc1?0)  at  several  wind  speeds.  Increasing  the  wind  causes 
the  sea  surface  to  resonate  more  strongly  around  the  peak  wavenumber  kp  while 
decreasing  the  peak  wavenumber  or  increasing  the  wavelength  of  the  peak  wave. 


logioW  *))  in  m  4) 


Figure  27:  Slices  of  Donelan-Pierson-Banner-Ryan  spectrum;  no  capillaries. 


The  sea  surface  admits  a  stochastic  Fourier  representation8: 


CM  =  » 


e+jkTr$(k)1/2G(k)dk 


m 


(15) 


where  (G(k)}  is  complex- valued,  zero-mean  Gaussian  white  noise  with  covariance  [1]: 


£[G(k)G(k')]  =  <5(k  —  k'). 


More  formally,  (G(k)}  is  a  tempered  random  held  that  is  filtered  by  <h(k)  to  gen¬ 
erate  a  Gaussian  sea  surface  [19].  Different  realizations  of  (G(k)}  generate  different 
realizations  of  the  sea  surface.  The  wavenumber  spectrum  <b(k)  is  based  on  Donelan 

8Allen,  J.  and  Frank.  J.  Ryan  [2000]  Stochastic  Sea-Surfaces  for  Radar  Simulations — Phase  1: 
Frozen  Gaussian  Surfaces,  preprint. 
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and  Pierson’s  [16]  wave  spectrum,  blended  with  Banner’s  [3]  observations  on  the  cap¬ 
illary  waves  and  the  spreading  function,  and  then  unified  by  Ryan9  for  both  a  fully 
developed  sea  and  a  young  sea. 

Figure  28  shows  a  digital  realization  of  a  Gaussian  sea  surface  plotted  over  a 
40  x40  nr  patch.  For  clarity,  the  vertical  scale  is  exaggerated.  By  a  digital  realization, 
we  mean  Equation  15  is  approximated  using  a  digital  filter  with  error  control.  For 
example,  the  wind  speed  is  set  as  U  =  10  (m/s).  Equation  14  predicts  a  significant 
wave  height  Hs  =  2.1376  (m).  Equation  13  estimates  Hs  =  2.1380  (m).  Enforcing 
this  close  agreement  between  the  theoretical  wave  height  and  estimated  wave  height 
is  one  example  of  the  error  control  and  suggests  that  the  sea-surface  simulation  has 
some  credible  features.  The  next  plot  shows  this  patch  in  a  larger  context. 


C(r):  *cut=30x  kf,  Nm= 1024 


Figure  28:  Gaussian  sea  surface. 


9Allcn,  J.  and  Frank.  J.  Ryan  [2000]  Stochastic  Sea-Surfaces  for  Radar  Simulations — Phase  1: 
Frozen  Gaussian  Surfaces,  preprint. 
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Figure  29  looks  down  upon  a  Gaussian  sea  over  a  2  x  2  km2  patch.  The  wind  is 
blowing  along  the  r\  axis  and  forcing  the  waves  to  line  up  approximately  orthogonal 
to  the  T\  axis.  If  a  transmitter  aims  along  the  r i  axis,  this  sea  surface  will  scatter  the 
bulk  of  the  rays  close  to  the  incident  plane.  This  scattering  is  the  in-plane  scattering 
of  Section  2.  This  in-plane  scattering  approximates  the  full-surface  scattering,  but 
needs  only  a  slice  C(ri?  0)  parallel  to  n  axis  for  scattering. 


Figure  29:  Top  view  of  Gaussian  sea  surface. 
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6  Scattering  From  Sea-Surface  Slices 


Figure  30  illustrates  the  scattering  model  for  the  sea-surface  slice.  The  incident  plane 
wave  reflects  off  the  flat  sea  surface  in  the  specular  direction  at  every  point.  As  the 
sea  surface  becomes  non- flat,  the  incident  plane  wave  reflects  off  the  local  tangent 
planes — subject  now  to  shadowing  and  blocking.  Shadowing  occurs  when  the  sea 
surface  prevents  the  incident  plane  wave  from  illuminating  a  reflecting  patch.  Block¬ 
ing  occurs  when  the  sea  surface  intercepts  a  reflected  ray.  Shadowing  and  blocking 
become  significant  at  small  grazing  angles. 


Figure  30:  (a)  Plane  wave  reflecting  from  flat  surface;  (b)  plane  wave  reflecting  from 
tangent  planes. 


35 


Figure  31  shows  a  slice  of  a  Gaussian  sea  surface  along  the  r i  axis.  An  antenna 
located  at  100  meters  has  a  1-meter  aperture  (green).  The  incident,  normal,  and 
reflected  rays  are  shown  for  those  reflectors  nearest  the  antenna.  This  plot  is  scaled 
so  that  the  angles  are  correct — the  angle  of  incidence  equals  the  angle  of  reflection 
measured  from  the  normal  to  the  surface.  The  rest  of  the  reflected  rays  come  from 
points  offscreen. 


Sea  Surface  Reflections:  9T=30 


range  (m):  U— 10  (m/s);  7^=  1  5C;  £1=0.83 


Figure  31:  Closeup  of  incident  and  reflecting  rays  that  arrive  in  the  antenna’s  aper¬ 
ture. 
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Figure  32  enlarges  the  view  to  show  a  long  slice  of  this  Gaussian  sea  surface. 
Plotted  on  the  sea  surface  are  scattering  points  (red  dots).  These  points  are  where 
the  transmitted  signal  reflects  off  the  sea  surface  and  arrives  in  the  antenna’s  aperture. 
That  is,  the  reflecting  patch  is  not  shadowed  from  the  transmitter  or  blocked  to  the 
receiver. 


o 

Sea  Surface  Reflections:  0T=3O 


Figure  32:  Gaussian  sea  slice  and  scattering  points. 

Associated  with  the  nth  scattering  point  are  the  following  reflection  parameters: 

•  Time  delay  A rn 

•  Elevation  angle  9n  from  the  antenna  pointing  to  the  scattering  point 

•  Grazing  angle  9gjn 

•  Reflection  coefficient  Tn  =  ro(0Sin) 
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These  scattering  points  model  the  delay-spread  function  by  summing  over  all  the 
reflectors: 

Nr 

h(r,r)  =  G](2(9t)6(t)  +  a(9T)  £  Gf(9n)TQ{9g,n)6(r  -  Arn) 

n=  1 

«  {of  (Or)  +  a(»T)  y  Gf  (0„)r„ J  «(r). 

The  diffuse  scatter  is  omitted  because  we  are  picking  up  this  effect  by  summing  over 
all  the  reflectors.  The  patch  factor  cl(9t)  blends  the  rough  surface  to  the  flat-surface 
model.  This  section  will  show  that  the  time  delays  are  negligible  for  a  narrow-band 
surface-level  antenna.  The  negligible  delays  allow  us  to  approximate  the  delay-spread 
function  like  multiplicative  noise.  In  this  case,  we  write  the  delay-spread  function  as 

Nr 

h(r)  «  GH2(9T)  +  a(9T)  £  Gf(9n)Tn  (16) 

n=  1 

to  act  on  the  signal  as  a  multiplication  operator, 

SR(t)  =  sT(t )  x  h(r(t)). 

Here,  r(t)  denotes  the  position  of  the  antenna  at  time  t.  Consequently,  the  multipath 
received  at  the  antenna  as  the  antenna  travels  across  the  sea  surface  at  various  speeds 
can  be  handled  by  resampling  h(r).  The  following  sections  compute  the  time  delays, 
angles,  and  the  patch  factor  a(9x)  that  blends  the  rough  surface  to  the  flat-surface 
model. 


6.1  Flat-Surface  Blending 


The  patch  factor  a(9x)  compensates  for  the  discretization  of  the  sea  surface  into  A r- 
sized  patches  and  blends  the  rough-surface  scattering  to  the  flat-surface  scattering  as 
the  sea  surface  gets  smoother.  Denote  the  upper  and  lower  heights  of  the  antenna 
aperture  as  hf,  and  respectively.  Let  A r  denote  the  increment  of  the  range 
samples.  The  simulations  all  use  the  plate  size  from  Didascalou  [15]: 


A  r  ~ 


A  t 

To' 


The  number  of  rays  that  could  reflect  off  the  “footprint”  of  the  antenna  looking  at  a 
flat  surface  is 


Nr 


round 


h-R  ~  hR  \ 
Artan($x)  j 


(U) 
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With  all  the  tiny  time  delays  set  to  zero,  the  flat-surface  multipath  scales  to  the  flat 
surface  discretized  with  A r  as 


Nr 


h(r )  =  G){2(0t)  +  N-1  ^i/2(-^T)r0,c(^T) 

71=1 

=  G'i/2(dT)  +  Gi/2(-dT)r0,c(^T). 

Consequently,  the  patch  factor, 

a(eT)  =  A 

nR  nR 

is  an  estimate  of  the  number  of  Ar-size  reflecting  patches  that  blend  rough-surface 
scattering  to  flat-surface  scattering  as  the  wave  heights  diminish. 


6.2  Reflecting  Points 

Figure  33  illustrates  the  geometry  required  for  a  ray  to  reflect  off  the  sea  surface  and 
arrive  in  the  antenna  aperture.  There  do  exist  back  reflections — reflections  that  come 
from  behind  the  antenna.  However,  if  the  transmitter’s  elevation  angle  is  small,  and 
the  sea  is  relatively  smooth,  back  reflections  into  the  surface-level  antenna  are  rare. 
Appendix  A  discusses  back  reflections  and  quantifies  this  claim. 


UN 


Figure  33:  Sea  surface  reflecting  a  ray. 


The  sea  surface  is  the  graph  of  £(r).  Points  on  the  sea  surface  have  the  form, 


Pc  = 


C(0 


Let  ujv  denote  the  upward  pointing  normal  at  p^.  The  unit  vector  pointing  to  the 
transmitter  is 

_  [  cos(dT) 

T  sin  (Ot) 
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The  angle  of  incidence  Oi  is  the  angle  between  u^r  and  u^: 


cos  (6i)  =  u^uT.  (18) 

By  adding  the  angles  in  Figure  33,  the  unit  vector  pointing  along  the  reflected  ray  is 


Ur  = 


cos  (0R) 
sin  (9r) 


(Or  —  Ot  +  2  Qi). 


For  the  reflected  ray  to  arrive  in  the  receiver’s  antenna  aperture, 


pr_  z 


(hR  <  hR  <  hR), 


there  must  be  a  dR  >  0  such  that 


Pr  =  Pc  +  dRuR 


or  that 


0  <  dR  = 


r  o  —  r 


^  ID  \  ■ 
cos(0R) 

Consequently,  the  reflected  ray  is  in  the  antenna’s  aperture  provided 

hR  <  C (r)  +  (r0  -  r )  tan(6»R)  <  h%.  (20) 

Equations  19  and  20  constitute  the  basic  reflection  algorithm.  The  only  modification 
is  that  if  a  ray  passes  these  tests,  it  still  must  be  tested  for  blockage  by  the  sea  surface. 


6.3  Time  Delay 

Figure  34  shows  the  time  delay  between  the  direct  and  reflected  rays  (and  is  distorted 
for  clarity).  The  reflected  ray  travels  the  diagonal  dR  from  the  point  of  reflection, 


CO")  ’ 


to  the  antenna 


Pi?  = 


This  reflected  ray  travels  the  distance 


dR  =  llPi?  —  Pell- 

At  the  same  time  the  reflected  ray  starts  for  the  antenna,  the  direct  ray  starts  from 
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Figure  34:  Time-delay  geometry. 


the  crossing  point  pA-  to  the  antenna  so  the  time  delay  is  the  difference, 


At 


dR  —  dT 


c 


The  crossing  point  is  located  at 


Pa-  =  Pc  +  dxel{6T+7r/2)  =  pfl  +  dTel0T. 


The  linear  system  is 

PC  -  P R 


so  that 


dTeidT  -  dxei{dT+n/2) 


cos  [6t) 
sin(0T) 


—  sin(0T) 
cos (8t) 


cos  (9t)  sin(#r) 

dj* 

sin(6lT)  —  cos  (6t) 

dx 

drp 


—  lx  det 


c 


-  r  o 


(r  -  r0)  cos(0T)  +  (C 


sin(0T) 

—  cos  (9t) 

-  hR)  sin(0r). 


The  time  delay  is 

At  =  c_1  |\/(r  -  r0)2  +  (C  -  hR)2 

~  {{r~  r0)  cos (6t)  +  (C  -  hR)  sin(0T))}  . 


(21) 


In  the  special  case  of  (  =  r0  =  0,  Equation  21  does  reduce  to  Equation  4.  How¬ 
ever,  Equation  21  is  an  excellent  example  of  the  numerical  difficulties  of  extracting 


41 


a  tiny  difference  by  subtracting  two  large  numbers.  The  following  implementation  of 
Equation  21  appears  to  overcome  this  numerical  problem — for  now.  For  brevity,  set 

A  r  =  r  —  r0,  Ah  —  £  —  Hr. 

Observe  that  Ah  is  negative  when  the  antenna  is  looking  down  on  the  reflector.  In 
this  notation, 

d,R  =  V  A  r2  +  Ah2. 

If  we  use  the  positive  angle  elevation  angle  from  the  antenna  to  the  reflector, 


-Ah 
A  r  ’ 

the  time  delay  can  take  the  forms: 


At 


c  1  {dR  —  (Ar  cos(9t)  +  A/isin(0T))} 
cos  ($T  +  0^)} 

(  A  Qc\ 


— {i  - 

c 

r,d'R  ■  2 

2 —  sm 


6.4  Reflected  Angles 


In  terms  of  Equation  16,  let  the  nth  reflector  be  located  at 

rn 

C (rn)  _ 

and  have  incident  angle  9i  n.  Because  the  incident  angle  is  measured  from  the  normal 
at  the  point  of  reflection  and  the  reflection  coefficients  are  functions  of  the  grazing 
angle,  the  reflection  coefficient  Tn  of  Equation  16  is  computed  as 


l-A  —  Tq  ,c(fig,n)i  dg,n  —  71"/2 


Likewise,  the  elevation  angle  0n  from  the  antenna  to  the  n  reflector  is 


9n  =  tan  1 


(  Q{rn)  -  hn\ 

V  rn  -  r0  )  ' 


42 


6.5  Scattering  Examples 


Figure  32  is  our  baseline  scattering  simulation.  Table  3  reports  on  the  scattering 
parameters  extracted  from  each  reflector.  We  use  these  parameters  to  model  the 
channel  in  Equation  16.  Because  the  antenna  is  fixed  at  a  single  position,  we  have 
only  a  “slice”  of  the  delay-spread  function. 


Table  3:  Reflection  data:  /  =  250  (MHz);  ay  =  0.5028  (m);  6t  =  30°. 


At  (ns) 

0n  (deg) 

|r0,c| 

7T 0,c  (deg) 

0g  (deg) 

8.768425 

-61.795469 

0.918591 

-3.989633 

45.798992 

8.897965 

-61.087893 

0.918344 

-4.000575 

45.217588 

10.387951 

-50.400428 

0.914882 

-4.154361 

39.138560 

11.387721 

-47.120155 

0.914587 

-4.167444 

38.740913 

11.539623 

-35.581093 

0.907388 

-4.486519 

31.739575 

11.677130 

-33.912428 

0.906298 

-4.534742 

30.967993 

11.973236 

-25.329341 

0.899442 

-4.837518 

27.041772 

11.855091 

-24.412361 

0.898195 

-4.892498 

26.456014 

13.256817 

-19.157287 

0.895087 

-5.029368 

25.121505 

12.497451 

-16.296921 

0.889205 

-5.287804 

22.983759 

12.617580 

-16.061290 

0.889020 

-5.295931 

22.923185 

13.620735 

-10.909591 

0.880003 

-5.690636 

20.350789 

14.290554 

-10.572835 

0.879618 

-5.707443 

20.254948 

13.660733 

-10.017306 

0.876356 

-5.849775 

19.479986 

14.924795 

-9.102366 

0.873784 

-5.961886 

18.911842 

16.353645 

-8.545893 

0.875865 

-5.871217 

19.368592 

19.623514 

-5.572819 

0.866038 

-6.298543 

17.392974 

21.318665 

-5.524093 

0.869955 

-6.128452 

18.128062 

48.720103 

-1.929041 

0.857431 

-6.671206 

15.973054 

100.354341 

-0.860771 

0.854128 

-6.813793 

15.488003 

155.560856 

-0.492800 

0.852544 

-6.882126 

15.265474 
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If  we  move  the  antenna  to  a  different  position,  a  different  channel  or  different 
“slice”  of  the  delay-spread  function  is  observed.  Figure  35  displays  all  the  channels 
obtained  by  moving  the  antenna  along  the  sea  surface.  Equivalently,  we  have  collected 
slices  of  the  delay-spread  function.  Consequently,  this  plot  shows  the  support  of  the 
“delay-spread”  function.  The  horizontal  axis  is  the  delay.  The  vertical  axis  is  the 
antenna’s  position.  As  the  antenna  moves  across  the  sea  surface,  the  nearby  waves 
offer  a  scintillating  effect,  while  the  more  distant  waves  offer  long-term  reflections 
delays  up  to  200  ns.  Even  these  longer  delays  are  not  resolvable  in  a  100-kHz  receiver. 


Sea  Surface  Delay-Spread:  0^=30 


Figure  35:  Range  and  delays  of  the  sea  surface;  transmitter’s  elevation  0t  =  30°. 
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For  comparison,  Figures  36,  37,  and  38  successively  decrease  the  transmitter’s 
elevation  angle.  The  large-scale  structure  of  the  sea  surface  is  lost  as  the  elevation 
angle  is  lowered — only  close-in  reflectors  are  operating. 


Sea  Surface  Delay-Spread:  0  =20 


100 


100  150  200 

delay  (ns):  t/=10(m/s);r  =15C;D=0.83 


Figure  36:  Range  and  delays  of  the  sea  surface;  transmitter’s  elevation  9t  =  20°. 
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Sea  Surface  Delay-Spread:  0r=lO 


Figure  37:  Range  and  delays  of  the  sea  surface;  transmitter’s  elevation  9t  =  10°. 
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Antenna  Position  (m) 


Sea  Surface  Delay-Spread:  0^=5 


Figure  38:  Range  and  delays  of  the  sea  surface;  transmitter’s  elevation  6 t  =  5°. 
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Table  4  lists  the  channel  from  the  6t  =  5°  simulation.  The  delays  are  not  in 
order  because  neighboring  reflectors  shoot  rays  that  land  anywhere  in  the  antenna 
aperture  covering  the  height  between  3  and  4  meters.  (See  the  closeup  in  Figure  31.) 
Consequently,  this  1-meter  aperture  “jumbles”  the  delays.  (The  same  effect  occurs 
in  Table  3.)  The  reflection  coefficients  are  still  relatively  large  because  the  grazing 
angle  is  between  2  to  10  degrees. 


Table  4:  Reflection  data:  /  =  250  (MHz);  cr^  =  0.5028  (m);  6t  =  5°. 


At  (ns) 

0n  (deg) 

To.c 

ZT0,c  (deg) 

0g  (deg) 

2.261470 

-13.246107 

0.776118 

-10.120476 

8.866531 

2.540322 

-10.765743 

0.770639 

-10.348577 

8.591683 

1.855261 

-9.942833 

0.732873 

-11.907600 

7.024713 

1.987390 

-9.211891 

0.728725 

-12.077468 

6.880722 

1.812893 

-7.967224 

0.702353 

-13.151415 

6.062877 

2.082713 

-6.274682 

0.695311 

-13.436511 

5.869202 

1.798164 

-3.639959 

0.636543 

-15.789676 

4.550052 

1.662531 

-2.565343 

0.591438 

-17.567117 

3.795084 

As  the  elevation  angle  decreases,  the  number  of  scatterers  decreases.  Although 
some  rays  were  still  arriving  when  the  elevation  angle  was  6t  =  2°,  dropping  the  ele¬ 
vation  angle  to  1  degree  eliminated  the  specular  reflection:  no  rays  were  reflected 
into  the  antenna  aperture.  At  this  low  grazing  angle,  only  the  wave  tops  are 
potential  scatterers.  This  restriction  of  the  scattering  to  the  wave  tops  significantly 
reduces  the  scattering  area.  Moreover,  the  sea  surface  is  discretized  at  increments  of 
size  A r  =  0.1  meter.  The  relatively  few  tangent  plates  that  approximate  the  ran¬ 
dom  wave  top  typically  will  not  precisely  align  in  both  height  and  angle  to  scatter 
into  the  antenna — provided  no  shadowing  or  blocking  occurs.  Therefore,  the  specular 
reflections  will  decrease  with  decreasing  elevation  angle. 
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7  Specular  Multipath  Noise  Simulations 


The  tangent-plane  scattering  model  is  applied  to  several  sea  surfaces  to  simulate  the 
multipath  noise  caused  by  specular  reflections.  Section  7.1  examines  the  multipath 
receive  by  the  Baseline  Quadrihlar  Antenna.  Section  7.2  applies  the  same  simulations 
to  the  OE-538  Low-Angle  Antenna.  The  simulations  show  the  following: 

•  The  specular  multipath  noise  migrates  to  the  upper  elevation  angles  as  the  sea 
becomes  rougher 

•  For  the  smoother  seas,  the  specular  multipath  noise  peaks  in  the  low  angles 

•  The  specular  multipath  rolls  off  at  very  low  grazing  angles 


7.1  Baseline  Quadrifilar  Antenna 


The  Baseline  Quadrihlar  Antenna  is  assumed  to  be  an  ideal  RCP  antenna.  In  the 
hat-surface  scattering  model  of  Equation  11,  this  antenna  receives  the  direct  path 
and  the  single  specular  reflection  as 

h(r)  =  Gr{9t )1^2  +  Gr{— ^t)1^2To,c(^t)- 


For  this  antenna,  the  direct-to-specular  ratio  (DSR)  from  this  hat  surface  is 


DSR(0t) 


Gr(6t) 

Gr(— 6it)|Toic(^t)|2 


Figure  39  plots  the  DSR  for  a  hat  sea  surface.  As  the  elevation  angle  Or  decreases 
to  zero,  Gr{±9t)  converges  to  Gr(0)  and  the  reflection  coefficient  vanishes.  Con¬ 
sequently,  the  DSR  becomes  unbounded  as  Or  decreases  to  zero.  This  same  effect 
appears  in  the  sea-surface  simulations,  but  is  confounded  by  shadowing  and  blocking. 

Figures  40,  41,  and  41  plot  the  specular  multipath  noise  scattered  from  the  sea 
surface.  The  multipath  noise  is  obtained  from  the  1-meter  aperture  of  Figure  32. 
The  antenna  pattern  of  Figure  13  resides  in  the  aperture.  The  antenna  starts  at 
r  =  100  and  moves  100  meters  in  0.1-meter  increments.  At  each  increment,  the 
multipath  noise  of  Equation  16  is  computed.  This  magnitude  of  the  noise  is  plotted 
as  a  function  of  range.  The  dB  plot  scales  the  noise  to  the  direct  path.  As  the 
elevation  angle  Or  decreases,  the  surface  appears  smoother  and  the  RCP  reflection 
coefficient  decreases. 
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Flat-surface  scattering  DSR:  hR=  3.5  (m) 


Figure  39:  Baseline  Quadrifilar  Antenna  DSR. 


Sea  Surface  Multipath:  9r=30  ',f= 250  MHz 


Figure  40:  Sea-surface  multipath:  transmitter  at  30°  elevation;  Baseline  Quadrifilar 
Antenna. 


50 


Sea  Surface  Multipath:  9r=10  ',f= 250  MHz 


Figure  41:  Sea-surface  multipath:  transmitter  at  10°  elevation;  Baseline  Quadrifilar 
Antenna. 


Sea  Surface  Multipath:  0r=5  \  f= 250  MHz 


Figure  42:  Sea-surface  multipath;  transmitter  at  5°  elevation;  Baseline  Quadrihlar 
Antenna. 
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Figure  43  displays  all  the  preceding  noise  plots  in  the  context  of  the  range  and 
elevation  angle  plane.  The  antenna  is  held  at  a  height  of  3  meters  and  moved  across 
the  sea  surface.  At  each  range,  the  elevation  angle  sweeps  over  0  to  30  degrees.  The 
noise  power — scaled  to  the  direct  path — is  plotted  as  a  function  of  range  and  elevation 
angle.  The  image  shows  a  relatively  strong  return  at  the  higher  elevation  angles  and 
little  return  at  the  lower  angles.  At  the  lower  elevation  angles,  the  grazing  angles  tend 
to  be  small.  Consequently,  the  sea  surface  reflects  less  of  the  incident  RCP  signal. 
Moreover,  at  very  low  angles,  almost  all  of  the  specular  reflections  are  blocked — the 
waves  either  shadow  the  scattering  patches  or,  after  reflection,  the  waves  block  the 
reflected  ray  to  the  antenna.  Consequently,  this  ray-tracing  model  has  the  specular 
noise  rolling  off  as  the  elevation  angle  is  decreased. 


Sea  Surface  Multipath  (dB  to  direct  path):/'  -250  MHz;  /?A=|3  4]  (m) 


100  110  120  130  140  150  160  170  180  190  200 

receiver  position  (m):  U=  10  (m/s);  T^=  15C;  £1=0.83 


Figure  43;  Sea-surface  multipath  for  the  Baseline  Quadrihlar  Antenna. 
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Figure  44  shows  that  when  the  wind  is  decreased  to  5  m/s,  the  sea  surface  becomes 
smoother  and  the  significant  wave  height  deceases.  The  large-scale  waves  still  line  up 
orthogonal  to  the  wind.  Consequently,  the  bulk  of  the  scattering  will  be  close  to  the 
incident  plane  and  is  approximated  by  in-plane  scattering  off  a  slice  of  the  sea  surface 
taken  along  the  n  axis. 


£(  r):  k  =30x  k 

J  cut  p 


-300  -200  -100  0  100  200  300 

r  (m):  U=5  m/s;  0=0°;  T  =15  C;  n=0.83;  H  =0.53439;  4c  =0.51699;  4  m‘/2=0.52357 

lv/  w  s  \  0 


Figure  44:  Sea-surface  realization  with  wind  at  5  m/s. 
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Figure  45  shows  the  corresponding  the  scattering  image  from  such  a  sea  surface — 
the  wind  is  only  moving  at  5  m/s  rather  than  the  10  m/s  of  Figure  43.  There  is  still 
blockage  and  shadowing  at  the  low  elevation  angles — but  the  effect  is  less  than  the 
rougher  sea  with  the  10  m/s  wind. 


Sea  Surface  Multipath  (dB  to  direct  path): 250  MHz;  /j^=[3  4]  (m) 


100  110  120  130  140  150  160  170  180  190  200 

receiver  position  (m):  U=  5  (m/s);  7/  =  15C;  £2=0.83 


Figure  45:  Sea-surface  multipath  for  the  Baseline  Quadrihlar  Antenna;  wind  at  5 
m/s. 
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As  the  wind  decreases,  the  sea  becomes  smoother.  Figure  46  shows  that  when 
wind  is  only  blowing  at  1  m/s,  the  scattering  is  nearly  uniform  in  range  and  follows 
the  DSR  curve  of  Figure  39. 
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-15 

-20 

-25 

-30 

100  110  120  130  140  150  160  170  180  190  200 

receiver  position  (m):  U=  1  (m/s);  r  ,=  15C;  12=0.83 

Figure  46:  Sea-surface  multipath;  Baseline  Quadrihlar  Antenna;  wind  at  1  m/s. 


Sea  Surface  Multipath  (dB  to  direct  path):  ,/=250  MHz;  «A=|3  4]  (m) 
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7.2  OE-538  Antenna 


The  flat-surface  scattering  model  of  Equation  12  consists  of  the  direct  path  and  the 
single  specular  reflection  that  splits  the  RCP  signal  into  RCP  and  LCP  rays  that 
arrive  at  a  non-ideal  antenna: 


h[r) 


C^rop^t)1/2  +  Gr,rcp(— 6it)1//2Po,rcp(^t) 

+G\r,lcp  ( — 9t)  1,/2ro,LCP  ($t)  • 


The  DSR  for  this  flat-surface  model  and  antenna  is 


DSR(6>t) 


_ Gr,b£?{6t) _ 

|Gr,rcp(— 9t  )1//2Po,rcp($t)  +  Gr,lcp(— ^t)1^2Po  ,lcp(0t)|2 


Figure  47  plots  the  flat-surface  DSR  for  a  sea  surface.  As  the  elevation  angle  Or 
decreases  to  zero,  the  reflection  coefficient  To.rcp^t)  converges  to  zero  so  that 


DSR(0°) 


_ (jr,rcp(0°) _ 

G^,lcp(0°)  |r0,LCp(0°)|2 


Flat-surface  scattering  DSR  for  the  OE-538:  hR=  3.5  (m) 


Figure  47:  DSR  for  the  OE-538 
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Ill  contrast  to  the  Baseline  Qnadrifilar  Antenna,  the  OE-538  is  still  responding 
to  the  cross-polarized  reflections  so  that  the  DSR  does  not  become  unbounded  as 
6t  decreases  to  zero.  This  same  effect  appears  in  the  sea-surface  simulations,  but  is 
confounded  by  shadowing  and  blocking.  Figures  48,  49,  and  50  show  the  specular 
multipath  for  the  wind  at  10,  5,  and  1  m/s.  Decreasing  the  wind  gives  a  smoother  sea 
surface.  The  smoother  sea  surfaces  migrate  the  multipath  noise  from  the  upper  ele¬ 
vations  angles  to  the  lower  elevation  angles.  The  smoothest  sea  reflects  the  multipath 
noise  essentially  as  the  DSR  curve  in  Figure  47. 


OE-538  Low  Angle  Sea  Surface  Multipath  (dB  to  direct  path):  f= 250  MHz;  hR=[ 3  4]  (m) 


100  110  120  130  140  150  160  170  180  190  200 


receiver  position  (m):  C/=  1 0  (m/s);  7A=15C;  0=0.83 
Figure  48:  Sea-surface  multipath  for  the  OE-538  Low-Angle  Antenna. 
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elevation  angle  (deg) 


OE-538  Low  Angle  Sea  Surface  Multipath  (dB  to  direct  path):  f= 250  MHz;  hR=[ 3  4]  (m) 


100  110  120  130  140  150  160  170  180  190  200 

receiver  position  (m):  U=  5  (m/s);  T  =15C;  £2=0.83 


Figure  49:  Sea-surface  multipath  for  the  OE-538  Low- Angle  Antenna;  wind  at  5  nt/s. 
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elevation  angle  (deg) 


OE-538  Low  Angle  Sea  Surface  Multipath  (dB  to  direct  path):  f= 250  MHz;  hR=[3  4]  (m) 
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100  110  120  130  140  150  160  170  180  190  200 

receiver  position  (m):  U=  1  (m/s);  r  ,=  15C;  =0.83 

Figure  50:  Sea-surface  multipath  for  the  OE-538  Low- Angle  Antenna;  wind  at  1  nt/s. 


59 


7.3  Multipath  Phase 


The  multipath  noise  is  at  baseband,  so  takes  on  complex  values.  The  preceding 
section  plotted  the  amplitude  as  a  function  of  the  receiving  antenna’s  range  and 
the  transmitter’s  elevation  angle.  This  section  shows  the  corresponding  plots  for  the 
multipath  phase.  Figure  51  shows  that  the  phase  is  relatively  stable  for  each  elevation 
angle.  The  exceptions  are  found  at  the  low  elevation  angles  where  the  magnitude  is 
small  and  no  phase  unwrapping  was  used  to  keep  the  angle  continuous.  Tables  3 
and  4  show  why  the  phase  is  relatively  stable — the  contribution  of  the  delays  to 
the  phase  is  negligible.  These  tables  also  show  that  the  grazing  angle  does  not  vary 
much — probably  because  of  the  narrow  aperture  of  the  receiving  antenna  and  its  short 
distance  from  the  sea  surface.  Because  the  grazing  angles  are  relatively  constant,  the 
phase  of  the  reflection  coefficient  is  also  relatively  constant. 


Baseline  QHA  Sea  Surface  Multipath  Phase  (deg):  f=25Q  MHz;  /zR=[3  4]  (m) 


100  110  120  130  140  150  160  170  180  190  200 

receiver  position  (m):  U=  10  (m/s);  7^=1 5C;  £2=0.83 
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Figure  51:  Multipath  phase;  baseline  QHA;  wind  at  10  m/s. 
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elevation  angle  (deg) 


Figures  52  and  53  plot  the  phase  of  the  multipath  as  the  wind  decreases  or  the 
sea  surface  gets  smoother.  These  plots  show  that  as  the  sea  surface  smoothes  out,  the 
phase  stabilizes  at  each  elevation  angle.  Consequently,  these  phase  plots  make  the 
obvious  observation  that  rougher  seas  cause  greater  phase  variations  than  a  smooth 
sea. 


Figure  52: 


Multipath  phase;  baseline  QHA;  wind  at  5  m/s. 
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elevation  angle  (deg) 


Baseline  QHA  Sea  Surface  Multipath  Phase  (deg):  f=25Q  MHz;  hR=[3  4]  (m) 
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Figure  53:  Multipath  phase;  baseline  QHA;  wind  at  1  m/s. 
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elevation  angle  (deg) 


7.4  Multipath  Correlation 


Figure  54  estimates  the  correlation  of  the  specular  multipath.  Each  elevation  angle  of 
the  transmitter  produces  a  multipath  noise.  The  plot  shows  the  correlation  between 
these  multipath  noises  as  a  function  of  elevation  angle.  Although  some  correlation 
is  visible  at  the  lower  angles,  there  is  also  little  power  at  these  angles.  In  the  upper 
angles,  the  multipath  is  decorrelated,  except  at  the  highest  angles. 

Baseline  QHA  Sea  Surface  Multipath  Phase  (deg) :  f= 250  MHz;  hR=[3  4]  (m) 
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Figure  54:  Multipath  correlation;  baseline  QHA;  wind  at  10  m/s. 
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elevation  angle  (deg) 


Figure  55  shows  a  decorrelation  in  the  low-elevation  angles  as  the  wind  decreases 
or  the  sea  surface  gets  smoother.  Correlation  in  the  upper  angles  increases. 


Baseline  QHA  Sea  Surface  Multipath  Phase  (deg):./' =250  MHz;  ^=[3  4]  (m) 
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elevation  angle  (deg):  U=5  (m/s);  7’r=15C;  £4=0.83 


Figure  55:  Multipath  correlation;  baseline  QHA;  wind  at  5  m/s. 
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Figure  56  shows  that  as  the  sea  smoothes  further  when  the  wind  decreases  to 
1  m/s,  the  multipaths  are  starting  to  correlate  at  the  upper  elevation  angles.  One 
conjecture  to  explain  these  correlation  plots  is  to  consider  a  single  small  plate  tan¬ 
gent  to  the  sea  surface.  Even  if  the  plate  reflects  into  the  antenna  aperture,  slight 
perturbations  of  the  plate  will  throw  the  reflecting  ray  out  of  the  antenna  aperture. 
Consider  now  a  collection  of  reflecting  plates.  On  a  rough  sea  surface,  this  set  can 
reflect  into  the  aperture  over  a  small  variation  in  the  transmitter’s  elevation  angle. 
A  large  shift  in  elevation  angle  illuminates  an  entire  different  set  of  reflectors — an 
uncorrelated  set.  Only  when  the  surface  is  sufficiently  smooth  do  large  patches  of  the 
sea  surface  become  correlated. 


Baseline  QHA  Sea  Surface  Multipath  Phase  (deg) :  f= 250  MHz;  /?A>=[  3  4]  (m) 
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elevation  angle  (deg):  U=  1  (m/s);  7/  =  15C;  £1=0.83 


Figure  56:  Multipath  correlation;  baseline  QHA;  wind  at  1  m/s. 
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7.5  Multipath  Envelope 


Identifying  the  envelope  of  noise  is  a  cottage  industry  for  the  communication  and 
radar  engineers.  This  envelope  identification  problem  forms  an  excellent  interface  be¬ 
tween  these  engineers  and  the  statisticians.  The  statistician’s  non-parametric  density 
estimators  require  an  immense  number  of  independent  samples.  The  envelopes  that 
the  communication  engineers  do  acquire  are  typically  narrow  band.  Consequently, 
the  independence  assumption  is  lost.  Therefore,  extending  the  non-parametric  den¬ 
sity  estimator  to  handle  narrow-band  time  series  is  an  excellent  research  topic. 

Figure  57  presents  an  Average  Shifted  Histogram  (ASH)  estimate  of  the  prob¬ 
ability  density  of  the  multipath  envelope  [39].  The  ASH  estimates  the  unknown 
probability  distribution  using  a  weighted  average  of  shifted  histograms.  The  plot 
shows  all  samples  of  {\h(r,  6*) | }  lumped  into  a  single  sample.  The  problem  with  this 
approach  is  that  each  “channel” — each  {\h(r,  0)|}  for  a  fixed  elevation  angle  9 — has 
a  different  power  level  that  confounds  the  identification  problem. 


Specular  Multipath  PDF: 4=250  MHz;  hR=[ 3  4]  (m) 


Figure  57:  Multipath  probability  density  estimate;  baseline  QHA;  wind  at  10  m/s. 
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Figure  58  is  an  estimate  of  the  root- mean-square  (RMS)  power  or  standard  de¬ 
viation  of  the  multipath  at  each  elevation  angle.  Power  is  lost  by  shadowing  and 
blocking  at  the  low  angles. 


Specular  Multipath  STD:  f= 250  MHz;  hR=\3  4]  (m) 


Figure  58:  Multipath  standard  deviation  estimate;  Baseline  QHA;  wind  at  10  m/s. 


To  reduce  the  effects  of  power  confounding  the  envelope  identification  problem, 
each  channel  is  normalized  by  dividing  by  its  estimated  standard  deviation.  Figure  59 
applies  a  Rician  goodness-of-ht  (GOF)  test  to  these  normalized  channels.  This  plot 
tests  if  these  normalized  channels  all  share  the  same  Rician  density  [36,  Eq.  5.59]: 


u 


p(u)  =  — a  exp 


crz 


u2  +  s2\  f  us\ 

) /o  W2)  ’ 


where  u  denotes  the  envelope,  s  denotes  the  magnitude  of  the  specular  component, 
and  cr2  is  the  variance  of  the  uncorrelated  noise.  These  parameters  are  estimated  from 
the  normalized  channels  and  appear  on  the  plot.  This  plot  also  provides  an  estimate 
of  the  specular-to-diffuse  ratio  or  K  factor  [36,  Eq.  5.60]: 

A'=101og10(Tj)  [dB], 
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Ricean  Goodness-of-Fit  Test:  P= 0% 


Figure  59:  Rician  GOF  of  the  normalized  multipath;  baseline  QHA;  wind  at  10  m/s. 

This  GOF  test  compares  a  histogram  against  the  Rician  density.  The  fit  is  good  for 
the  large  magnitudes  but  deviates  from  Rician  in  the  small  magnitudes.  The  GOF  is 
the  observed  significance  level  P  that  measures  how  closely  the  histogram  follows  the 
Rician  density.  Small  values  of  the  observed  significance  level  (P  <  0.1)  reject  the  fit. 
However,  a  key  assumption  for  the  significance  test  is  that  the  samples  be  independent. 
Dropping  the  independence  assumption  is  a  basic  research  problem  for  GOF  theory. 
No  data  pruning  or  whitening  was  applied  to  this  multipath.  Consequently,  the  low 
significance  is  confounded  by  the  correlations  in  the  multipath. 

Because  of  the  good  fit  over  the  large  envelope  values,  it  is  natural  to  remove 
the  small  magnitudes  from  the  multipath  and  test  again.  Specifically,  the  multipath 
channels  with  9  <  5°  were  omitted  from  the  normalized  multipath.  Figure  60  shows 
a  GOF  applied  to  this  multipath  with  the  low  angles  removed.  The  small  deviations 
from  Rician  are  omitted,  but  the  fit  is  still  not  tight  on  the  Rician  density. 

Similar  results  can  be  obtained  with  the  Weibull  and  the  //-distributions.  This 
blind  fitting  is  a  pitfall  of  non-parametric  identification.  That  is,  it  is  very 
tempting  to  throw  densities  at  the  multipath  and  hope  to  “get  lucky”  with  a  large 
significance  level.  In  contrast,  some  clutter  research  mixes  physics  and  statistics  to 
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Ricean  Goodness-of-Fit  Test:  P= 0% 


Figure  60:  Rician  GOF  of  the  normalized  multipath — angles  below  5°  omitted;  base¬ 
line  QHA;  wind  at  10  m/s. 


arrive  at  new  GOF  test  and  distributions  [12].  A  fully  rational  approach  asks: 

Given  the  Gaussian  sea  surface  and  specular  scattering,  what  envelope 
density  will  we  see  based  on  physical  reasons? 

This  question  is  part  of  the  larger  question  that  asks  how  the  sea-surface  statistics  (i.e., 
the  variance  of  the  wave  heights,  the  surface  correlation  length,  the  slope  statistics) 
map  to  the  multipath  statistics. 
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8  Modeling  Issues 


This  model  provides  qualitative  simulations  of  specular  sea-surface  multipath  imping¬ 
ing  upon  a  surface-level  antenna.  This  ray  model  assumes  the  following: 

•  Ray  propagation 

•  Narrow-band  receiver 

•  Transmitter  and  receiver  aligned  with  the  wind 

•  Fully  developed  sea 

•  Signal  wavelength  approximately  1  meter 

•  Local  tangent-plane  scattering 

•  No  surface  aerosol  effects 

•  No  atmospheric  effects 

•  Far-held  antenna  pattern 

The  surface-level  antenna  and  narrow-band  receiver  force  the  time  delays  of  each  ray 
into  the  nanosecond  regime.  A  receiver  with  a  kilohertz  bandwidth  cannot  register 
these  time  delays.  Consequently,  the  time  delays  can  be  ignored. 

Aligning  the  transmitter  and  receiver  with  the  wind  over  a  fully  developed  sea 
puts  the  rays  orthogonal  to  the  bulk  of  the  waves.  Consequently,  most  scattering  is 
in-plane.  In-plane  scattering  lets  us  scatter  from  a  “slice”  of  the  sea  surface  rather 
than  scattering  over  the  entire  surface. 

The  transmitter’s  wavelength  A t  is  approximately  1  meter.  Consequently,  scat¬ 
tering  from  sea-surface  features  that  are  less  than  0.1  meters  in  size  can  be  ignored. 
Because  the  capillary  waves  are  millimeter  size,  the  rough-surface  coefficient  is  approx¬ 
imately  1  and  limits  the  diffuse  scattering.  For  this  reason,  specular  reflecting  patches 
are  limited  to  0.1  meters  in  size.  The  tangent-plane  scattering  does  not  model  the 
curvature  effects  of  the  surface  nor  the  reflections  off  the  top  of  the  waves.  Likewise, 
the  effects  of  the  surface  aerosols,  whitecaps,  bubbles,  and  spray  are  omitted. 

The  atmosphere  is  assumed  to  deliver  plane  waves  to  the  scattering  region  around 
the  antenna.  Consequently,  atmospheric  scintillation  and  ducting  are  omitted.  Only 
straight-line  ray  propagation  is  used. 

Finally,  antenna-to-surface  coupling  effects  are  ignored  by  keeping  the  antenna 
at  least  a  wavelength  away  from  the  local  surface.  These  simplifications  allow  a 
straight-forward  model  description: 
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Each  tangent  plate  on  the  surface  m,ay  reflect  a  specular  ray  into  the  an¬ 
tenna  aperture,  provided  the  ray  is  not  shadowed  or  blocked. 

This  simple  model  provides  a  qualitative  description  of  the  specular  multipath 
effects  for  a  surface-level  antenna.  In  light  of  these  approximations,  it  makes  little 
sense  to  polish  a  small  model  part  in  the  presence  of  the  following  large-scale  effects 
and  approximations.  First,  the  antenna  pattern  changes  as  the  antenna  gets  closer 
to  the  sea  surface.  The  simulations  of  this  report  kept  the  base  of  the  antenna  3 
meters  above  the  mean  sea  surface.  The  largest  significant  wave  height  was  2  meters, 
which  allows  the  waves  to  come  within  a  wavelength  of  the  antenna.  The  sea  surface 
will  modify  the  antenna  pattern.  For  example,  a  tilting  ground  plane  will  distort  the 
main  beam.  Consequently,  the  near-held  effects  of  the  sea  surface  must  be  computed 
whenever  the  antenna  is  within  a  wavelength  of  the  sea  surface. 

Second,  two-dimensional  propagation  over  the  sea  surface  “slice”  needs  to  be 
compared  to  3-D  propagation  over  the  full  sea  surface.  However,  comparing  ray 
models  does  not  address  the  quality  of  the  simulations.  The  real  comparison  should 
be  between  the  ray  models  and  high-quality,  partial-differential  equation  models. 
For  example,  the  Variable  Terrain  Radio  Parabolic  Equation  (VTRPE)  code  handles 
the  rough-surface  scattering  over  a  slice  of  a  sea  surface  [37].  Comparing  the  two- 
dimensional  ray  model  against  the  VTRPE  simulation  provides  a  practical  assessment 
of  the  ray  model  propagating  over  a  slice.  If  this  comparison  is  “good  enough”  for 
communication  simulations,  extending  the  ray  model  to  3-D  propagation  is  justified. 

Third,  this  multipath  simulation  is  essentially  local.  We  assume  that  plane  waves 
carry  the  signal  into  the  scattering  region.  This  plane-wave  assumption  implies  that 
there  is  no  distortion  by  ducting  or  ionospheric  scintillation — even  at  low  elevation 
angles.  This  assumption  is  questionable  and  the  effect  on  the  signal  needs  to  be 
quantified  and  compared  to  the  multipath  effect.  For  example,  measurements  of  low- 
elevation  SATCOM  signals  over  a  sea  surface  have  been  carried  out  at  the  Scripps 
Pier  in  La  Jolla  by  Ken  Anderson,  SSC  San  Diego  [2],  Figure  61  reports  on  the 
received  SNR  as  a  function  of  increasing  satellite  range  or,  equivalently,  decreasing 
elevation  angle  as  the  GPS  satellite  heads  down  toward  the  horizon.  Reexamination 
of  these  measurements  or  conducting  a  similar  held  experiment  provides  a  low-cost 
measurements  that  can 

•  quantify  the  plane-wave  assumption 

•  recover  the  multiplicative  noise 

•  benchmark  the  multipath  simulations 
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Scripps  Pier  GPS  Observation 


Normalized  Subsatellite  Range  (km) 

097631622  Fit  Dec.  24,  1998 

Figure  61:  Measured  GPS  satellite  SNR  at  a  sea-level  receiver. 

Fourth,  the  antenna  measurements  need  to  be  upgraded  to  include  the  effective 
length  or  effective  height  h  of  the  antenna  [32,  Chapter  14].  The  effective  length 
determines  the  received  voltage  Vq  at  the  antenna’s  terminals  as 

Ro  =  hTE4, 

where  Ej  is  the  incident  electric  held  plane  wave  [40].  If  the  electric  held  is  composed 
several  plane  waves,  say 

Ej  =  E^i  +  Ej,2  +  Ej,3, 

the  effective  length  the  antenna’s  response  to  each  plane  with  the  correct  relative 
phase  is 

Vo  =  hTEfil  +  hTE,2  +  hTEj)3. 

Therefore,  getting  the  effective  length  by  the  Numerical  Electromagnetic  Code  or 
measurements  is  necessary  for  multipath  simulations. 

Finally,  regardless  of  how  the  multipath  is  computed,  there  is  fascinating  research 
linking  the  sea-surface  statistics — the  standard  deviation  of  the  height,  the  correlation 
length,  the  distribution  of  the  slopes — to  the  statistics  of  the  specular  multipath. 
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A  Back  Reflections 


These  simulations  use  only  reflections  from  the  sea  surface  in  front  of  the  surface-level 
antenna.  Here,  “front”  refers  to  the  region  between  the  receiver  and  the  transmitter. 
If  the  transmitter’s  elevation  angle  is  small  and  the  sea  is  relatively  smooth,  almost 
no  rays  ever  reflect  from  behind  the  antenna.  This  appendix  quantifies  this  claim  so 
that  these  rare  events  can  be  omitted. 

Figures  62  and  63  display  reflections  arriving  from  behind  the  antenna.  Both 
figures  show  a  relatively  steep  wave  face  is  needed  to  bounce  a  low-angle  ray  back 
into  the  surface-level  antenna.  The  relevant  question  is — how  common  are  such  steep 
waves? 


Figure  62:  Back  reflection — upper  ray. 


Recall  that  £(r)  denotes  the  sea  surface  and  £(?"!,  0)  denotes  a  slice  of  the  sea 
surface.  The  slopes  of  the  waves, 

C'0t,°)  :=  J^Ot,0), 

have  a  Gaussian  distribution.  For  the  sea  surface  with  the  wind  of  10  m/s,  (/(r^O) 
is  zero  mean  with  a  estimated  standard  deviation 

<j£/  ~  0.1616. 

Assume  a  back  reflection  shown  in  Figure  62  occurs  at  sea  level  with  a  transmitter  at 
30  degrees  elevation.  The  wave  slope  required  typically  exceeds  —75  degrees  but  the 
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Figure  63:  Back  reflection — lower  ray. 


smallest  wave  slope  is  —38  degrees.  The  probability  of  such  a  slope  is  small: 

Prob[tan(slope)  <  — 40°]  =  Prob[slope  <  atan(— 38°)]  =  2.1  x  1CT5. 

The  probability  that  no  back  reflections  of  this  type  happen  over  the  100  meters 
behind  the  receiver’s  antenna  exceeds  98% — assuming  that  no  shadowing  or  blocking 
occurs.  Consequently,  if  the  transmitter’s  elevation  angle  is  small  and  the  sea  is 
relatively  smooth,  the  back  reflections  are  ignored. 
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